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Abstract 

fS| ■ The formalism to describe the scattering of a weakly bound projectile nucleus by 

a heavy target is investigated, using the Uncorrelated Scattering Approximation. 
The main assumption involved is to neglect the correlation between the fragments 
of the projectile in the region where the interaction with the target is important. It 
is shown that the angular momentum of each fragment with respect to the target is 
conserved. Moreover, when suitable approximations are assumed, the kinetic energy 
of each fragment is also shown to be conserved. The 5-matrix for the scattering 
of the composite system can be written as a combination of terms, each one being 
^T* ! proportional to the product of the 5-matrices of the fragments. 

'. 

Key words: Nuclear Reactions, Scattering Theory, Three-Body Problem, Halo 
Nuclei, Elastic Scattering, Inelastic Scattering, Breakup Reactions. 
PACS: 24.10.Eq, 24.50. +g, 03.65.Nk, 25.10.+S, 25.70.Bc, 25.70.Mn 

X' 



(N 



1 Introduction 



In recent years there has been much interest in the description of halo nuclei. 
These are weakly bound and spatially extended nuclear systems, so that some 
fragments of the nucleus (generally neutrons) have a high probability of being 
at distances larger than the typical nuclear radii (see Refs. [1,2] for a general 
review on these nuclei). 
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Reactions induced by halo nuclei on different targets are a useful tool to inves- 
tigate the structure of these nuclei. A proper description of the reaction mech- 
anism of the collision of a halo nucleus with a target is a complex quantum 
mechanical problem, which requires the explicit treatment of the continuum 
of breakup states. This treatment can be performed by discretizing the con- 
tinuum of breakup states. However, the price to be paid is that the relation 
between the scattering observables measured and the structural properties 
involved in the problem becomes troublesome. 

When the scattering energy of the projectile is sufficiently high, the interaction 
of the halo nucleus with the target can be considered to occur in such a short 
time that the fragments of projectile are practically "frozen" in given relative 
positions. This "frozen halo", also known as "adiabatic" approximation, is the 
basis of different approximate treatments which have been successfully applied 
to the interpretation of halo nuclei scattering measurements [3,4,5]. 

The validity of the "frozen halo" approach depends on the comparison of the 
collision time Tc and the internal time 77. The collision time can be estimated 
as Tc = a/v, where a is a typical length scale of the interaction and v is the 
relative velocity. The internal time is given by Tj = h/e where e is a typical 
excitation energy of the projectile. Following this argument, the "frozen halo" 
approximation will be valid provided that the collision time is much shorter 
than the internal time tq <C 77. 

Nevertheless, tidal forces which tend to distort the structure of the projectile 
while interacting with the target, are also very important. This distortion has 
a time scale td which can be estimated in terms of the tidal force F T = Vr/a 
through the expression a = l/2F T /mr|,. Thus, the validity of the "frozen halo" 
approximation requires also tq <C td- 

In this paper we focus on a different approach named the Uncorrelated Scat- 
tering Approximation (USA), that we first introduced in [6]. Its basic approxi- 
mation consists in assuming that the interaction between the fragments of the 
halo nuclei can be neglected during the collision time. Thus, the fragments 
of the projectile scatter separately from the target, each of them with a frac- 
tion of the total energy of the projectile. The validity of the USA is linked 
to a correlation time t r = K/A, where A measures the energy spread of the 
ground state of the projectile when the interaction between the fragments is 
neglected. So, the USA is expected to be valid when tq <C tr, but, in contrast 
with the "frozen halo" approach, it does not depend on the strength of the 
tidal forces. 

The adiabatic or "frozen halo" approximation implies that the relative coordi- 
nates within the projectile are constant during the scattering process. As we 
shall show, the USA implies that the relative angular momentum between each 
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fragment of the projectile and the target is conserved. Moreover, the kinetic 
energy of each fragment is also conserved. 

In section 2 we develop the formalism of the uncorrected scattering approx- 
imation (USA) making use of the expansion of the three-body T-matrix. In 
section 3 we present an application of the USA approach to the case of elastic 
scattering and breakup of deuteron on 58 Ni. In section 4 the conclusions are 
presented. 



2 The Uncorrelated Scattering Approximation from a T-matrix 
approach 



We consider a two-body system (A, B) which is scattered from a heavy target 
(T). Neglecting spin and other internal degrees of freedom of the fragments 
and the target, the Hamiltonian is given by 



„ P 2 f 

H = — + - + v AB + v AT + v BT 

p2 p2 

= TT^ - + + Vab + Vat + Vbt ' W 

2m A 2,m b 

where M = m A + m B and \i = m A m B j{m A +m B ). The interactions v AB , v A t 
and v B t depend on the relative coordinate and momenta of the interacting 
particles. 

To define a scattering process, we consider an asymptotic situation in which 
the fragments are very far from the target. In this case, the asymptotic Hamil- 
tonian is given by 

p2 p2 

^=2M + Y, +VAB =2M + h - < 2 > 



The eigenstates of H corresponding to a given total energy E which are rele- 
vant for a scattering process can be expressed as a free wave \ki, LM L ) on the 
relative coordinate times an eigenstate I Mr) of the internal Hamiltonian h. 
Note that these states can be characterized by a given angular momentum and 
projection. It should be noticed that the eigenstate \4>i, I Mi) corresponding to 
an eigenvalue of h, may be a bound or a continuum state. The momentum ki 
is related to the energy of the state through % 2 kf /2M + e,i = E. We can couple 
the product states to a total angular momentum J, to give kf, L, J, JMj). 
Since spin projection Mj is conserved, it will be dropped in the following 
derivation. The effect of the interaction v A t + v B t of the projectile with the 
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target is given by a T-matrix, connecting the eigenstates of H with energy 
E. Formally, the T-matrix can be expressed as 

T(E) = vat + v B t + {vat + v B t) 1 TT T(E) , (3) 

where -E 4 " = E + ie. 

The scattering amplitude is proportional to the on-shell matrix elements con- 
necting the initial bound state \<f> g , k g ; L, I, J) with a final state \<f>f, kf] V , I', J), 
which may be bound or unbound. Here, the quantities k g and kf represent the 
relative momenta of the two-body system (A, B) with respect to the target T 
in the initial and final channels, respectively. The ^-matrix can be expressed 
as 



{(f) f , k f ; L', I', J\S\<f>g, k g ; L, I, J) = 

= 6 f , g + 2m( ( f> f , k f ; L', /', J\T{E)\4> g , k g , L, I, J) . (4) 

Then, in order to apply the Uncorrelated Scattering Approximation, we make 
use of the fact that the interaction vat + vbt conserves separately the angular 
momenta La and Lb of the fragments of the projectile with respect to the 
target. Note that the T-matrix will not, in general, conserve La and Lb, 
due to the interaction vab between the fragments, which is included in the 
propagator. However, a reasonable approach is to neglect the term vab in the 
propagator, while we rescale it by a factor A, which can be a function of J, and 
will be determined later. This approximation implies dismissing the relatively 
weak interaction vab in the region where Vat + Vbt is strong (it contributes 
to second order to the T-matrix). Hence, we make the approximation 

1 ~ A (K) 

E + ie-H E + ie-Ho' K ' 

where H = Ka+Kb is the Hamiltonian which contains only the kinetic energy 
terms. The factor A included in this expression may be complex, describing 
in that case the loss of flux due to the excitation of the target and/or the 
fragments of the projectile. 

The above approximation leads to the expression 



T{E)^f{E) (6) 
f(E) = \{v AT + v BT ) + X(vat + vbt) f+ 1 gv f{E) . (7) 

Note that T(E) is the T-matrix for the scattering of two uncorrelated particles 
A and B with a renormalized interaction \{vat + vbt)- 
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The eigenstates \<pi, kf, L, I, J) of H Q can be expanded in terms of eigenstates 
of#o, 

\4>ii h] L, I, J) — J dagi(a)\a,Ei(a);L,I,J). (8) 

In the above expressions, Ei(a) = (E — Cj)sec 2 (a) is the eigenvalue of H 
with a the hyperangle which depends on the ratio of the internal momentum 
q and the relative momentum hi, i.e., q(ki,a) = ki^J/i/M tana. The explicit 
expression of gi(a) is given in the appendix A. 

The matrix elements of T(E) become 



^ f ,k fl L',l',J\f(E)\<P 9 ,k 9 ,L,I,J) 
= J jda'dag f {a')g g {a){a\E f (a')-L'j\J\f(E)\a,E g (a)-LJ,J). (9) 

Thus, Eq. (9) reflects that the on-shell matrix element of T, between eigen- 
states of H corresponding to the energy E, are given in terms of off-shell 
matrix elements of T between eigenstates of Ho corresponding to different 
energies E g (a), Ef(a'). 

Besides the assumptions involved in (9), in what follows we make a further 
approximation by replacing the off-shell matrix elements of T(E) by on-shell 
matrix elements. So, while we retain the a dependence of the matrix elements, 
we substitute the energy dependence for the on-shell value. Thus, the matrix 
elements in (9) can be approximated by 



(a', E f (a'); L', I', J\f(E)\a, E g (a); L, I, J) 
~ \A f (a')(a', E; L', I', J\T(E)\a, E; L, /, J)A 9 (a) . (10) 

The factors A f (a') and A 9 (a) are introduced to ensure unitarity. Note that the 
typical values of E^a) — E are of the order of the kinetic energy of the relative 
motion of the fragments, which is related to the correlation time through U/tr. 
The operator T(E) connects states in a range of energies which is determined 
by the collision time through AE = H/tc- So, this on-shell approximation is 
justified provided that the correlation time is long compared to the collision 
time. 

Now, for the purpose of evaluating the on-shell matrix elements of T(E) we 
can make use of an expansion in hyperspherical harmonics, which allows us 
to write down the operator T(E) between states with definite values of the 
angular momenta La and Lb of each fragment with respect to the target. Note 
that these magnitudes are conserved by the interaction. After some algebra 
(see Appendix B for details), we finally get 
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(<f> f ,k f -L'j',J\T{E)\<t> g ,k g -LJ,J) 
-EE (<l>f\K')vi>(K\<l> a )u(L A ,L B \LJ) K j(L',r\L A ,L B ) K ,j 

KK> L A ,L B 

x (K', E; L A , L B , J\T(E)\K, E; L A , L B , J) (11) 

with 

(K', E; L A , L B , J\f(E)\K, E; L A , L B , J) 
= J J d/M/7/£ >Lfl (/?)^^ 

(12) 

where the functions fL A L B (fl) an d the overlaps (K\<j) g )Li are given in the 
appendix B, and {L A , L B \L, I) KJ are the Raynal-Revay coefficients [7]. The 
hyperangular variable (3 determines the partition of the energy E between 
the fragments A and B of the projectile, so that e A = Ecos 2 (f3) and e B = 
Esm 2 (P). 

Expanding the T-matrix up to third order and using the results of appendix 
C the matrix elements of the operator T(E) can be written as: 

((?, E\f(E)\(3, E) = 5((3 - (3'){t A (e A ) + t B (e B ) + 2mt A {e A )t B {e B )} , (13) 

where t A (e A ) denotes a on-shell two-body matrix element for the A — T scat- 
tering (analogously for t B ). 

Recovering the angular momenta, the matrix elements of the T-matrix be- 
tween states with given K, L A , L B , J values, result 

(K', E; L A , L B , J\f(E)\K, E; L A , L B , J) = J d(3f^ LB ((3)f^ LB ((3) 
x ^t A (E cos 2 (3) + t B (E sin 2 f3) + 2mt A (E cos 2 /3)t B (E sin 2 (3)^, (14) 

which, according to Eq. (4), can be written in terms of S, i.e., the S'-matrix 
which describes the scattering with boundary conditions given by H , as 

(K', E; L A , L B , J\S\K, E; L A , L B , J) 
= J d(3f« A , LB {(3)f^ LB {(3)S A {E cos 2 (3)S B {Esin 2 (3) . (15) 

If the values of K, K' are restricted by a maximum hyperangular momentum 
K m , the above integral can be approximated by quadratures. It is particularly 
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convenient to consider the quadrature associated to the Jacobi polynomials 
that determine the functions fL A L B (P) with a number of points N equal to 
the number of K values allowed. Explicitly: 



(K',E;L A ,L B ,J\S\K,E;L A ,L B ,J) 

N 

= E( K »L A ,L B (n\K) LA , LB S A (e n A )S B (e B ), (16) 



n=l 



where, following [6,8], the coefficients (K\n) l a ,l b are given by 



< A »— = J^W^W^ fLL ° iM ■ (17) 

The quadrature points (3 n correspond to the zeros of the function fL A ,L B (ft)- 
The explicit expression for the quadrature weights w n can be found in [6]. Note 
that, in terms of the quadrature points, the individual energies of particles A 
and B are given by e A = Ecos 2 ([3 n ) and e B = Esin 2 (f3 n ). 

Then, collecting these results, we may write 



(<l> f ,k f ;L',I',J\S{E)\<l> a ,k g ;L,I, J) 
- H ((pf,kf]L\l',J\n]L A ,L B ,J)S A (e r X)S B (e B )(n]L A ,LB,J\(pg,k g ]L,I, 

L A ,L B ,n 



with 

((pi, h; L, I, J\n\ L A , L B , J) = ^{(t>i\K) LJ {L, I\L A , L B ) K j(K\n) La ,l b - 

K 

(19) 

To finish with this analysis we proceed by evaluating the parameter A which 
determines the renormalization of the interactions v A and v B . For this purpose, 
we write the individual ^-matrices in terms of the phase shifts 

S A (e n A ) = exp (2i5(\v A ; L A , e^j ; S B (e n B ) = exp ^2i8(Xv B ; L B , e n B )j . 

(20) 

In many approximations the phase shift is found to be proportional to the 
potential. Under this situation, the phase shift scales with the parameter A. 
This is for instance the case of the eikonal and WKB approximations, where 
the phase shift is obtained as an integral of the potential along a classical 
trajectory. In general, there is not a simple relation between the phase shift 
and the potential, but only for the purpose of estimating in a simple way 
the value of A, we will assume that this relation holds approximately. So, 
5(\v A ; L A , e A ) = \5(v A ; L A , e A ). 
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In the limit of weak interactions, one expects that breakup effects are negli- 
gible, and thus the elastic S matrix should be well described by the folding 
interaction. Let us assume that the ground state has / = 0, implying L = J. 
Then, the folding interaction, which depends on the relative coordinate and 
momentum is given as 

v F = (<fr g \v AT + V BT \(j) g ) (21) 

and the expression for the elastic ^-matrix due to the folding interaction can 
be written as 

S F {E, J) = exp (2i5{v F ; J, E)j . (22) 
Thus, comparing Eqs. (22) and (18) in the limit of weak interactions, one gets 

5(v F ; J, E) = XS L ij , (23) 

where 

hu= X! \(n;L A ,L B , J\(j) g , k g ;L, I, J)\ 2 l5(v A ; L A ,e A ) + S(v B ;L B ,e B ) > ■ 

n,L A ,L B 

(24) 

Using this result, the product of S'-matrices can be written as 

S A (e n A )S B (e n B ) = exp ^7rp(n, L A , L B , J)5(v F ; J, E)^j , (25) 

where 

p(n, L A , L B , J) = ^5(v A ; L A , e n A ) + 5(v B ; L B , e n B ) }/6 L u- (26) 

Note that p(n, L A , L B , J) is a real, positive number, with typical values around 
unity. It fulfils the relation 

p(n, L A , L B , J) | (n; L A , L B , J\(f> g , k g ; L,I,J)\ 2 = 

= ]T \{ n] L A ,L B ,J\<P g ,k g -LJ,J)\ 2 = l . (27) 

The value of p(n, L A , L B , J) measures how strong is the interaction of the 
fragments with the target in the configuration which is characterized by the 
quantum numbers n,L A ,L B , compared with the average interaction over all 
the configurations that contribute to the elastic channel. Note that Eq. (25) 
implicitly states that any absorption effect, due to the excitation of the target 
or the fragments of the projectile, that can be described by the use of a complex 
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folding potential, will also appear in the uncorrelated three-body calculation, 
scaled by the factor p(n, L A , Lb, J). 

We would like to stress that Eq. (18) can be applied to elastic, inelastic and 
breakup scattering. In particular, the elastic S'-matrix elements are given by 



{<j> g ,k g ;L,I,J\S(E)\(j>g,kg;L,I,J) = 
]T |(n;L A ,L B ,J|0„A; s ;L,/,J)| 2 ^(e^)^( e S). (28) 

La, Lb,™ 

The inelastic ^-matrix to an excited bound state \<p e ) is 



{<t> e ,k e ;L'J', J\S(E)\<f> g ,k g ;L,I, J) = 

( n ~, l a, L B , J\<t> 9 , k g ; L, I, J) (0 e , k e ; L', I', J\n\ L A , L B , J)S A (e A )S B (e 

La, Lb,™ 



It is also possible to calculate the ^-matrix elements for the breakup leading to 
specific states of the continuum, characterized by an asymptotic linear internal 
momentum p and quantum numbers {L',I',J}. According to Eq. (18) these 
are given by 



S J l I;L 'Ap) = (<P(p), k{p); L', I', J\S{E)\4> g , k g ; L, /, J) 
= (<f>(p),k(p);L',I',J\n;L A ,L B ,J) 

LA,L B ,n 

x (n; L A , L B , J\<j> g , k g ; L, I, J)S A {e n A )S B {e n B ) . 



(30) 



Note that, due to energy conservation, h 2 k 2 (p)/2M + h 2 p 2 /2fx = E. The dif- 
ferential breakup cross section is calculated directly from the above ^-matrix 
elements by means of the expression 



d^iiMi'ip) 7T2J+1 



dp 



k 2 g 2I+l 



(31) 



An appealing feature of our previous formulation of the USA, presented in 
[6], is that it leads to compact expressions for the integrated breakup cross 
sections. Although the derivation presented in this work differs to some extent 
from our previous formulation, it also leads to similar closed formulae for 
the breakup cross sections. In particular, the total integrated breakup cross 
section to all the final continuum states characterized by the angular momenta 
L', V, J is given by (c.f. with Eq. (36) of Ref. [6]): 
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K„ 



o*(LI^L'n = £—±-{ £ \(K,L>,I>,J\S\(i> 9 ,k g ;m,J)f 

K g Zl -t- 1 L K=L , +I , 

- E I (<l>e, K\ L', I', J\S\(j)g, kg] L, I, J) | 2 

e 

- <5l,z/<M (<f> g , k g ; L, I, J\S\4> g , k a ; L, I, J)| 2 }. (32) 

This expression is readily obtained by summing upon V and V in Eq. (31) 
and applying the closure relation (B.9) to the final states. 

It is also possible to obtain a compact expression for the breakup cross section 
corresponding to a total angular momentum J, <7j\ This is achieved upon 
summation of a b j u (LI — > VI') on the angular momenta V and I' and taking 
into account the completeness property of the states \K, V , I' , JMj) . This 
leads to the closed expression 



Tli = ^Jf^j{(^ k ^ L > J > J \S + S\<P 9 , k g ; L, I, J) 



- ]T \((f) e ,k e ;V,I',J\S\(f) g ,k g ;L,I,J)\ 

e,L>,I> 



2 



-\{<Pg,kg;LJ,J\S\<P g ,k g ;LJ,J)\ 2 Y (33) 

Then, within the USA, the integrated breakup cross section for a given total 
angular momentum is calculated as the dispersion of the operator S in the 
ground state of the projectile, subtracting the contribution of the other bound 
states. 



In summary, in this section we have derived an approximate formula to eval- 
uate the elastic, inelastic and breakup ^-matrix elements corresponding to 
a three-body scattering process. These approximated S'-matrix elements are 
expressed in compact form by Eq. (18). The simplicity of this formula re- 
lies on the fact that it relates the complicated three-body collision matrix 
to a simple superposition of two-body ^-matrices, weighted by some factors 
which depend on analytical coefficients and the structure of the composite. 
Expression (18) provides also a possible physical interpretation for the states 
\n; La, L b , J). These states are eigenstates of the Hamiltonian H in a dis- 
crete basis. According to (18), the initial state is decomposed in the full set of 
states \n; La, Lb, J)- Upon neglection of the inter-cluster interaction vab, the 
two particles evolve in this set of states, interacting with the target through 
the interactions va and vb, but preserving the individual energies and angular 
momenta of the constituents. The superposition of states resulting after the 
interaction with the target is finally projected onto the physical final state for 
the process under consideration, denoted V, I', J). 
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3 Application to deuteron scattering on Ni 

As an illustrative example, in this section we apply the USA method to the 
reaction d+ 58 Ni at £^=80 MeV. This reaction has been previously analyzed 
by several authors within the CDCC scheme [9,10,11,12], providing a satisfac- 
tory agreement with the existing elastic data [13]. Therefore, in this work 
we have adopted the CDCC calculation as the benchmark result to com- 
pare our method with. In order to obtain the elastic and breakup observ- 
ables, we performed CDCC calculations similar to those reported in [9,10]. 
The proton-target and neutron-target interactions were described in terms 
of the Becchetti-Greenless global parameterization [14]. The p-n interaction, 
which is required to construct the deuteron ground state and continuum bins, 
was parametrized as v(r) = t> exp(— r 2 /a 2 ), with t> = —72.15 MeV and 
a = 1.484 fm [9]. Following the standard procedure, the p-n continuum was 
discretized into energy bins. The partial waves 1=0, 2, 4, 6 and 8 were included 
in the CDCC model space. The odd partial waves were not considered, since 
their influence on the dynamics turned out to be negligible. The maximum 
excitation energies for these waves were e max =40 MeV (s- waves), 45 MeV (d 
and /-waves), 50 MeV (g-waves) and 60 MeV (h- waves). For simplicity, the 
proton and neutron spins are ignored. These CDCC calculations were carried 
out with the computer code FRESCO [15]. 

The USA calculations were performed using Eq. (18). According to this expres- 
sion, the three-body ^-matrix is constructed by superposition of the on-shell 
cluster-target ^-matrices, evaluated at the appropriate energies and angular 
momenta. For computational convenience, we found useful to evaluate these 
S'-matrices using the WKB approximation. This has the extra advantage that 
it makes exact the scaling property of the phase shift with the parameter A, 
thus making more sensible the estimate (23) for this parameter. 



3.1 Elastic scattering 

Since the USA method is formulated in terms of the ^-matrix, we first com- 
pare the elastic S'-matrix elements in both approaches. The ground state wave- 
function was described with the same Gaussian potential used for the CDCC 
calculation. 

In Fig. 1 we present the CDCC and USA calculations in the form of an Ar- 
gand plot, in which the real and imaginary parts of the ^-matrix elements 
are represented in the x and y axis, respectively. Note that, since the internal 
spin is neglected, the orbital angular momentum for the projectile-target rela- 
tive motion L coincides with the total angular momentum J. The CDCC and 
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Re (S L ) 

Fig. 1. Argand plot for the CDCC (filled circles), USA (open circles) and clus- 
ter-folded (crosses) calculations, corresponding to the d+ 58 Ni elastic scattering at 
E d =80 MeV. 

USA calculations are represented by filled and open circles, respectively. For 
comparison purposes, we have also included the cluster-folded (CF) calcula- 
tion (crosses), in which the projectile-target interaction is given by the sum 
of the proton-target and neutron-target interactions, folded with the deuteron 
density, according to Eq. (21). Note that the latter is equivalent to a CDCC 
calculation without continuum bins. Therefore, the difference between the fold- 
ing and the CDCC calculations evidences the effect of the breakup channels 
on the elastic scattering. From the curves in Fig. 1, it becomes apparent that 
the presence of the continuum reduces the modulus of the ^-matrix for all the 
partial waves, due to the loss of flux to the breakup channels. It can be seen 
that the USA method also succeeds on reproducing this deviation from the 
CF calculation, giving a result close to the CDCC. 

The differential elastic angular distributions derived from the above ^-matrices 
are shown in Fig. 2, along with the data of Stephenson et al. [13]. The dotted 
line represents the folding calculation, which clearly overestimates the data, 
particularly at intermediate angles. Inclusion of the continuum within the 
CDCC scheme (dashed line), produces a significant reduction of the elastic 
cross section and improves significantly the agreement with the data. Notice 
that, despite this improvement, the data are still somewhat overestimated at 
angles around 60°. Two different USA calculations are presented in Fig. 2. The 
first one, represented in the graph by the dotted-dashed line, was performed 
with a maximum hyperangular momentum K m =30. This calculation agrees 
fairly well with the data in the whole angular range. Note that for c m . < 25° 
the agreement with the data is excellent. At intermediate angles, the USA 
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Fig. 2. Elastic scattering angular distribution for the d+ Ni reaction at 
Ed=80 MeV. The meaning of the curves are indicated by the labels. 

reproduces well the interference pattern, although the absolute normalization 
is somewhat underestimated. This suggests that the method will overestimate 
the breakup cross sections, as it will be discussed below. For cm. angles above 
80° the USA calculation shows some oscillations, which are not present on nei- 
ther the CDCC calculation nor the experimental data. These spurious wiggles 
are partially suppressed when the cutoff hyperangular momentum is increased 
to K m =40, as shown by the solid line. We found that for higher values of K m 
the angular distribution remains basically unchanged. Thus, we can conclude 
that for K m =40 we have convergence of the USA calculation. 



3.2 Breakup reaction 



We next analyse the breakup channel for the same reaction. Due to the absence 
of experimental data, we will compare our model directly with CDCC. 

We first study the integrated breakup cross section, as a function of the to- 
tal angular momentum J. In the USA calculation, we can make use of the 
closed expression (33). The results are shown in Fig. 3. It can be seen that 
the USA calculation (solid line) reproduces the qualitative behaviour of the 
CDCC (dashed-line). In particular, both distributions predict a maximum of 
the breakup cross section at J = 17. However, our model yields significantly 
more breakup than the CDCC, as we anticipated from the behaviour of the 
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Fig. 3. Integrated breakup cross section, as a function of the total angular momen- 
tum, J for the reaction d+ 58 Ni at 80 MeV. The dashed and solid line correspond to 
the CDCC and USA calculations, respectively. 

elastic angular distribution. 

As described in the previous section, the method here developed provides 
also compact expressions for the partial integrated cross sections leading to 
continuum states with definite values of {L',I',J}, according to Eq. (32). Using 
the same physical ingredients as in the elastic scattering, we have applied 
this formula to calculate the integrated breakup cross section for I'=0 12. 
The results are displayed by open circles in Fig. 4. It becomes clear that 
the main contribution to the breakup cross sections comes from the partial 
waves I'=0 and I'=2. Also included in this figure are the CDCC results (filled 
circles). Note that within CDCC these integrated cross sections are obtained 
as a result of large set of coupled equations, in which all possible values of 
V and I' need to be coupled simultaneously. From this figure it becomes 
apparent that the discrepancy of the breakup cross section between the two 
methods is mainly due to the components J'=0 and I'=2. For I' > 2 the 
two methods are in excellent agreement. This result becomes more notable if 
one compares the simplicity of Eq. (32) with the large set of coupled-channels 
equations involved in the CDCC. Furthermore, when high excitation energies 
and angular momenta are to be included, convergence problems typically arise 
in CDCC. This is in contrast with Eq. (32), which allows one to calculate the 
breakup to highly excited states and large values of I' in a simple way without 
the numerical difficulties involved in a CDCC calculation. 

Besides the integrated breakup cross sections, we can also calculate the breakup 
to specific states of the continuum. The ^-matrix element describing the tran- 
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Fig. 4. Integrated breakup cross section, as a function of the final p-n angular mo- 
mentum, I'. The filled circles (joined with dashed lines to guide the eye) correspond 
to the converged CDCC calculation (I max = 8), while the open circles joined by solid 
lines are the USA prediction obtained with the closed expression (32). 



sition from the ground state to a continuum state with asymptotic linear mo- 
mentum p and angular momenta {L 1 , 1', J} is given by Eq. (30). Two different 
models where used to describe the continuum states within the USA method. 
In the first one, we used a discretized continuum in terms of continuum bins, 
as in the CDCC calculation. Since these bins are normalizable the overlaps 
{(pf\K) LI appearing in Eq. (30) are calculated in exactly the same way as for 
the ground state. In a second calculation we used for the final states the true 
scattering wavefunctions, calculated at a given excitation energy. Note that 
these states are no longer normalizable and so special care has to be taken in 
order to calculate their overlaps with the hyperspherical basis. Details of the 
evaluation of these overlaps can be found in Appendix A. In this calculation, 
we found convenient, although not essential, to work with analytic expressions 
for the wavefunctions. In particular, for the ground state we used the Hulthen 
wavefunction: 



with a = 0.2317 fm" 1 and (3 = 7a. For the continuum states, we used the so- 
lution of a separable potential, whose ground state coincides with the Hulthen 
wave function, thus guaranteeing the orthogonality with the ground state 
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wavefunction. These are explicitly given by 




1 



e ipr + f(p) 



e ~ipr _ e ~/3r 



} 



(35) 



(2tt) 3 



r 



with 



/(p) = - /3- 



/? 2 + p 2 (/3 2 + p 2 ) 2 
2p 2p(a + PY 



(36) 



where a and /3 are the same as for the ground state. These wavefunctions are 
normalized as {4> p >\4> p ) = S(p' — p). Note also that this potential acts only on 
the s-waves; for 1^0, the continuum sates are simply given by planes waves. 

In Fig. 5, we compare the modulus of the breakup ^-matrix for a total angular 
momentum J = 17 which, according to Fig. 3, corresponds to the maximum 
of the breakup distribution. The continuum states with I' = and I' — 2 
have been considered for the comparison. In the latter, the separated contri- 
butions for L'=15, 17 and 19 are shown. Two different CDCC calculations 
are presented. The first one, represented by open circles, uses a model space 
with I = 0, 2 only. However, this model space is not enough to achieve con- 
vergence of the ^-matrix elements. In analogy with [10], we had to include 
partial waves up to / m ax=8. The results of this CDCC calculation in the aug- 
mented space is shown by the filled circles in Fig. 5. The USA calculations 
with the Hulthen potential are represented in Fig. 5 by the solid lines. Finally, 
the USA calculation performed with the continuum bins is represented with a 
histogram, to emphasize the fact that this calculation uses a discretized con- 
tinuum. Both calculations used a maximum hyperangular momentum K m =30. 
With K m =AQ the results are only slightly changed. In general, we find a good 
global agreement between USA and CDCC. For small values of p (i.e. low ex- 
citation energies) the two methods yield very similar results. On the contrary, 
for large values of p all our calculations tend to predict more breakup than 
the CDCC. Interestingly, the USA distributions exhibit a bump at p ~1 fm -1 
(i.e. e ~ 41 MeV), which is not observed in the converged CDCC with / max =8, 
but appears nevertheless in the CDCC with J max =2. The major disagreement 
between the USA and CDCC occurs for the breakup to s states. In this case, 
the two methods give very different predictions for p > 0.3 fm -1 . 

The differential breakup cross section is readily calculated from the above S- 
matrix elements using Eq. (31), and summing upon V and J. The results are 
displayed in Fig. 6. The filled circles correspond to the CDCC calculation in 
the model space with 7 ma x=8. The solid line is the USA calculation with the 
p-n separable interaction. The histogram corresponds to the USA calculation 
in the discretized continuum, generated with the Gaussian potential. We see 
that, at small excitation energies, the USA and CDCC calculations are in very 
good agreement. At higher excitation energies, the USA predicts systemati- 
cally more breakup, as expected from the behaviour of the ^-matrices analysed 
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Fig. 5. Breakup S'-matrix elements for a total angular momentum J = 17 and fi- 
nal p-n relative angular momentum I'=0 and 2. The filled circles correspond to 
the CDCC calculation with s and d waves only, whereas the open circles are the 
CDCC calculations with the model space 7 = 0,2, 4, 6, 8. The solid line is the USA 
calculation with analytic wavefunctions derived from a separable potential. The his- 
togram correspond to the USA calculation using the bins constructed with Gaussian 
potential. 

above. For V = 6, the USA calculations obtained with the continuum bins 
agree nicely with the CDCC. In the cases J'=0 and I'=2 we show also the 
adiabatic coupled-channels calculation by Amakawa and Tamura [16]. This 
calculation follows a similar trend to that of the CDCC, which is understood 
by recalling that the adiabatic approximation can be regarded to an approxi- 
mated CDCC calculation in which the internal Hamiltonian is replaced by a 
constant [9]. 

From the analysis performed in this section, we conclude that, despite its 
formal simplicity, the model developed in this work provides a good description 
of the reaction observables. In particular, we have shown that the method 
describes fairly well the elastic and breakup scattering of d+ 58 Ni at £^=80 
MeV. Our comparison with the CDCC calculations suggests that the USA 
tends to overestimate breakup to highly excited states on the continuum. 

In comparing the CDCC calculation with the USA results, it should be taken 
into account that the USA approach is not simply an approximation to the 
CDCC calculation. The latter is performed assuming that the interactions 
between fragments and target, that are complex, can be approximated by 
local potentials. This determines the off-shell behaviour of the interaction 
and, in particular, prevents coupling to highly excited states of the projectile. 
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Fig. 6. Breakup excitation function for the reaction d+ 58 Ni at Ed=80 MeV. The filled 
circles represent the CDCC calculations. The solid lines are the USA calculations 
with the Hulthen wavefunctions. The histogram is the USA calculation obtained 
with a Gaussian potential and a discretized continuum. The dotted-dashed lines are 
the adiabatic calculation of Amakawa and Tamura [16]. 

In contrast, the USA approximation relies on the fact that off-shell matrix 
elements of the interaction are substituted by on-shell ones. This seems to be 
the reason behind the larger breakup cross sections to highly excited states. 



4 Summary and conclusions 



In this paper we have revisited the uncorrelated scattering approximation 
(USA) originally introduced in Ref. [6]. We reviewed the basic assumptions 
involved within the USA approach and analyzed its capability to describe 
elastic and breakup scattering reactions. In what follows we summarize the 
main ingredients and results obtained in this work. 

The description of the scattering reaction mechanism provided by the USA 
model is based on three basic approximations. First, in the case of a weakly 
bound system interacting with a heavy target, we ignore the interaction be- 
tween the fragments in the three-body T-matrix propagator. Thus, the opera- 
tor T(E) may be replaced by T(E)/X, being T(E) the uncorrelated T-matrix 
for a renormalized interaction. This has the property to conserve the angu- 
lar momentum of each fragment of the projectile with respect to the target. 
However, it results essential to conserve the interaction in the asymptotic 
states, which are eigenstates of the Hamiltonian H . Expanding these states 
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into eigenstates of Hq (Hamiltonian containing only kinetic energy terms), one 
finally ends up with the evaluation of the off-shell T(E) matrix elements. 

The second approximation involved is to replace the off-shell matrix elements 
of T(E) by the on-shell ones. This is justified because the range of off-shellness 
in the matrix element is small compared with the energy range of the operator. 
However, a direct substitution of the off shell matrix element by on-shell ones 
may lead to the breakdown of the unitarity property. Thus, in this work we 
have made use of a renormalization operator, which is formally equivalent to 
the Democratic Mapping procedure described in [17]. The matrix elements of 
T(E) between eigenstates of H can be evaluated using an expansion in hyper- 
spherical harmonics. This allows us to transform analytically the asymptotic 
states into states with definite values of the angular momenta La and Lb of 
each fragment with respect to the target, and take advantage of the fact that 
these magnitudes are conserved by the interaction. 

The third approximation is to expand up to third order the three-body op- 
erator T(E) in terms of the three-body operators Ta{E) and Tb{E) which 
contain only the interaction of one of the particles with the target. Then the 
on-shell matrix elements of T(E) can be evaluated, and the result shows that 
the operator T(E) does not connect states in which the total energy E is dis- 
tributed differently between particles A and B. Moreover, the matrix element 
of T(E) is completely determined by the on-shell matrix elements of the two- 
body operators t^e^) and ts(efi) evaluated at the corresponding energies of 
the two particles. 

Our analysis shows that the ^-matrix describing the scattering of a composite 
system by a target, in the uncorrelated scattering approximation, is given as 
a combination of products of the S'-matrices describing the scattering of the 
fragments evaluated at the corresponding energies and angular momenta. The 
renormalization factor A can be obtained for each J value to ensure that, in 
the weak coupling limit, the USA reproduces the elastic scattering calculated 
with the folding interaction. This renormalization factor allows us to include 
the effect of excitation of the target and/or the fragments of the projectile, 
which is essential in nuclear collisions, in a way which is fully consistent with 
a complex folding interaction. 

The USA keeps also some resemblance with other impulse approximations 
recently applied to the scattering of weakly bound nuclei [18,19]. In analogy 
with the multiple scattering of the T-matrix method presented in [19], we 
start with an approximated expansion of the few-body T-matrix, in which 
the inter-cluster interaction is neglected in the propagator. However, while 
in [19] the derivation is performed in a linear momentum representation, we 
choose a partial wave description. Moreover, although both methods express 
the total T-matrix in terms of two-body amplitudes, the approximations that 
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lead to the two-body T-matrices are quite different. 

There exists also some formal analogy between our main result, Eq. (18), 
and the semiclassical Glauber approximation [3], in the sense that in both 
methods the scattering amplitude depends on a superposition of the product 
of the individual ^-matrices of the fragments. However, despite this apparent 
similitude, we stress that the approximations involved in both approaches 
are very different. First, unlike the Glauber model, the USA does not make 
any semiclassical assumption. Furthermore, the Glauber model is based on 
the frozen halo or adiabatic approximation (i.e., it neglects the excitation 
energies of the internal Hamiltonian), while the USA neglects the inter-cluster 
potential (the so called impulse approximation), but retains the full kinetic 
energy operator. 

Finally, it should be stressed that the main differences found between the USA 
and the CDCC calculations arise from the fact that the USA calculation gives 
rise to larger breakup cross sections to states with high excitation energies. 
This result might be related to limitations of the USA approach, but it could 
also be due to the fact that the CDCC calculations assume local interactions 
between the fragments of the projectile and the target, and this determines 
the off-shell nature of the interactions. Accurate experimental measurements 
of breakup cross sections at high excitation energies would surely help to draw 
more definite conclusions on the validity of the local, momentum independent, 
complex interactions used in the CDCC approach or, by contrast, on the 
reliability of the presence of relevant off-shell components in the interaction 
as suggested by the USA calculations. 



A Expansion of the channel wavefunctions in terms of the hyper- 



The bound states of the projectile, as well as the normalizable bins of contin- 
uum states, can be expressed in momentum representation as 



where the state \q\IMj) is normalized so that {q';IMj\q;IMj) = 5(q — q'). 
The corresponding channel states can be written as 



angle 




(A.l) 




(A.2) 
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It is convenient to write this state in terms of states with given hyperangle a 
and kinetic energy Ei(a), which are given by 

which leads to 



\(f>i, h; L; I , J) = J dagi(a)\a,E(i,a);L;I,J) (A. 4) 

g l {a) = <PM{h,a))l^^j 1 ' 2 ■ (A-5) 
The last term arises from the normalization (a', E\a, E) = 5(a' — a). 
The continuum states with asymptotic momentum p are given by 

\<j>[p), fc(p); /, L, J) = b, fc(p); /, L,J)+ J dqg^-)(p, q)\q, k(p); I, L, J) , (A.6) 

where the first term in the RHS represents the plane wave and the second term 
represents incoming waves. The integral has a pole at q — p, with a residue 
R(p) and a principal part, so 

100), Hp); I, L, J) = [1 + mR{p)\ \p, kip); /, L, J) 

+ V I dqg(-)(p,q)\q,k(py,I,L,J). (A.7) 

This state can be written in terms of the hyperangle as 



\<P(p),k(py,I,L,J) = [l + t7rR(p)](^-) ' \a p ,E;I,L,J) 



da P/ 



+ V J dag { -)(p,q(a)) (^j \a, E(p, a); I, L, J) , 

(A.8) 



where 



p [m h 2 P 2 h 2 k( P ) 2 /4 x 

tan (a- = J — , E = — — + - ;y , A.9 

v PJ k(p) V /i 2fi 2M ' y ' 

q [M _ , h 2 q 2 h 2 k(p) 2 lk x 
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B On-shell matrix elements in the hyperspherical basis 



In this appendix we show in detail the procedure used to evaluate the on-shell 
matrix elements of the operator T(E). First, we make use of an expansion in 
hyperspherical harmonics which allows us to write down 

/ dag{(a)A\a)\a, E; L, I, J) = £<tf|&) LJ |tf, E; L, /, J) , (B.l) 

J K 

where we introduce the states 

\K, E; L, I, J) = fj 2 daf^(a)\a, E; L, I, J). (B.2) 

The functions fL,i( a ) are given by 

/&(«) = ^( CO s«) L+1 (sina) /+1 P^&|(cos2a), (B.3) 

where P n 2 ' 2 are the Jacobi polynomials of degree n and are some 
normalization constants, whose explicit expressions can be found in [6]. 

The coefficients {K\(pi)u are explicitly given by 



(KMli = E ^m^K', K) (B.4) 

K' 

d LI (K') = J daf%(a)g!( a ) (B.5) 
X LI {K\ K) = J dafZ(a)A\a)f^(a) . (B.6) 

To get these expressions we have made use of the closure property of the 
hyperspherical harmonics 

E =5(a-a>). (B.7) 

K 

Then, we end up with the following expression for the T-matrix elements 



(<j> f ,k f ;L',r,J\T(E)\<f> g ,k g ;L,I,J)~ 

E (<f>f\K') VI , (K\(j> g ) LI (K', E; L', I', J\f(E)\K, E; L, I, J) . (B.8) 

K',K 

The requirement that the approximated expression preserves unitarity for her- 
mitian interactions leads to 

Y,( K '\4>i)Ll(4>i\K)LI = &K,K>. (B.9) 
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On the other hand, if the interaction is constant, the matrix elements of T(E) 
will be diagonal and independent on K. In this case, the T-matrix should not 
couple different internal states. This leads to 

EW,M*/I*)l' = <W- (B.10) 

K 

These conditions can be achieved by an adequate choice of A l LI (K',K). In 
what follows we shall assume, for definiteness, that in each spin channel / 
there is at most one bound state, while the rest of the states correspond to 
the continuum. So, for the bound state in channel / we take 

A b LI (K,K')=6 K , K/ /^/p[n (B.11) 

with P^j = J2k I c l/(-^)| 2 ) while for the rest of the states in channel / we 
consider a unique symmetric matrix A C LI (K',K) which fulfils 

^A c LI (K,K')c b LI (K') = 0, 

K' 

£ A C LI (K, KMIAK', K[) £ cUK,)d LI {K[) 

K U K[ i+b 



Note that this procedure is equivalent to orthogonalize all the continuum states 
with respect to the ground state, and then apply the "democratic mapping" 
procedure to the continuum states. 

It should be noticed that the sum in K is extended in principle up to infinity. 
Note that in this case P\ T — 1. In practice, the sum is taken up to a maximum 
value K m , which is obtained to get convergence in the calculations. So, in 
general, P\ jl < 1, specially for the higher partial waves, for which L is close 
to K m . Note that, in contrast to what was done in Ref. [6,20], no physical 
meaning is attached to the parameter K m , which should be taken as large as 
possible, until convergence is achieved. 

In order to make use of the fact that T(E) conserves the angular momenta of 
the fragments, one can use the Raynal-Revai transformation [7]: 

\K,E;L,I,J) = £ \K,E;L A ,L B ,J)(L A ,L B \L,I) KJ . (B.14) 

La,L b 

The states with different K values can be expanded to get states with definite 
values of the kinetic energy of each fragment. We can write, in terms of the 



(B.12) 

c b u(K)c b LI (K>) 

'K,K> ^ • 

(B.13) 
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hyperangle, 

\K,E;L A ,L B ,J) = JdPf* AiLB (P)\P,E;L A ,L B ,J). (B.15) 

This state can be characterized in terms of a product state of particles A and 
B with energies e A = Ecos 2 ({3) and e B = Esin 2 ((3). 



C Multiple scattering expansion of the T operator 

In this appendix we evaluate the on-shell matrix elements (/?', E\T(E)\(3, E), 
where the state \(3, E) is an eigenstate of H = K A + K B , corresponding 
to energies e A = E cos 2 (3 and e B = Esin 2 (3. Analogously, the state \/3',E) 
corresponds to e' A = E cos 2 (3' and e' B = Esin 2 (3'. Note that E = e A + e B = 
e' A + e' B . We can write explicitly 

\(3,E) = N(f3)\e A ,e B ), (C.l) 

where N{/3) = \j E sin(2^) is the square root of the Jacobian of the transfor- 
mation from {e A ,e B } to {{3,E}. 

We consider the multiple scattering expansion of the operator T(E) up to 
third order 

T(E)=T A (E)+T B (E) 

+ T A (E)G (E)T B {E) +T B {E)G (E)T A (E) 

+ T A (E)Go(^)T B (^)Go(^)T A (E)+T B (E)Go( J B)T A ( J B)G'o( J B)T i3 (E), 

(C.2) 

where Gq(E) = (E + — K A — Ks)^ 1 is the free 3-body propagator and 

T A (E) = Xv AT + Xv AT G (E)T A (E) . (C.3) 

The matrix operator T A (E) can change the kinetic energy of particle A, 
through its off-shell components, but it can not modify the kinetic energy 
of B. Hence, we can write 

(P',E\T A (E)\P,E) =S([3-[3')(e A \t A (e A )\e A ), (C.4) 

where we have introduced the operator 

tA(e A ) = Af at + \v AT g A {e A )t A (E), (C.5) 

with g A {e) = (e + — K A ) the free 2-body propagator. Note that the operators 
T A and t^ differ on the kinetic operator appearing in the propagator. While T A 
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is defined with the full kinetic energy operator, i.e., K A + Kb, the propagator 
in tA contains only the kinetic energy operator associated with particle A, 
Ka- Therefore, Ta should be understood as a three-body operator, whereas 
tA corresponds to a two-body operator. For simplicity, in the following, we 
drop the energy argument of the three-body operators Ta, T b , Gq, which is in 
all cases E, but we retain it in the two-body operators. 

The contribution of the second order terms can be expressed as the sum of a 
pole term and a principal part 

((?, E\T A G T B \P, E) = m5{{3 - (3')t A (e A )t B (e B ) 

+ (ef A \tA(e' A )\e A )(e' B \t B (e B )\e B ) 
E - e' B - e A 
E\T B G T A \P, E) = mS(/3 - (5')t A {e A )t B {e B ) 

+ <e^|tB(e^)|6B>(e^|t A (e A )|e^) 
E — e A — es 

where we have introduced the short notation t A (e A ) and ts(es) for the on-shell 
matrix elements (e^lt^e^le^) and (eB|ts(es)|es), respectively. 

As for the third order terms we have 



((3',E\TaG TbG T a \(3,E) 

= f de' A (e' A \t A (e' A )\e' A )—-^ 7f r (e' B \t B (E - e" A )\e B ) 

J bj^ - e A - e B 

x F+ \ - &\tA(eA)\e A )N(P)N(p) . (C.8) 
t,^ — e A — e B 

To evaluate this contribution we make use of the following identities 

g A (e')g A (e) = (^(e') - g A (ej) (C.9) 

t A (e') - t A (e) = t A (e') (g A {e') - g A {e)}t A {e) . (CIO) 

The second of these expressions indicates that the operator t B (E — e" A ) can be 
approximated by t B (e B ), and the difference would be of higher order in the 
T-matrix expansion. Thus, one may write 

(P',E\T a G T b G Ta\P,E) 

~ (e B |t B (e B )|e s )(e^|t A (e^)^(e^)^(e A )t A (e A )|e^)7V(/3)7V(/3') . (C.ll) 



Making use of the above identities, one gets 
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(I3',E\T A G„T B G„T A \S,E} 

e A -e A 

A similar derivation for the other third order term gives 

(/3',E\T B G a T A G T B \P,E) 

^ (e' A \t A (e A )\e A ) ^Mfk>lf£> - ^ B (e B )\e B ) , 

&b - e' B 

So, the sum of the two third order contributions reduces to 

((3',E\T A G T B G T A \(3,E) + ((3' , E\T B G T A G T B \(3, E) 
= Wb |tg(e B ) \e B ) (e' A \t A (e' A )\e A ) N ^ N ^ 
e A — e' A 

+ (e^|t A (e A )|e A )(e^|t i? (e- B )|e i? ) iv( ^ )Ar(/3/) ^ 

which cancels exactly the principal value part of the second order terms, 
Eqs. (C.6) and (C.7). Therefore, collecting these results, we have that the 
on-shell matrix elements of the T operator up to third order can be written 

as 

((?, E\f(E)\(3, E) = 5((3 - (3'){t A (e A ) + t B (e B ) + 2mt A {e A )t B {e B )} . 

(C.15) 
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